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Galaxies in the Fields of 2; ~ 1.5 Radio-Loud Quasars 
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Patrick S. Osmer,^ Donald W. McCarthy,^ Aaron S. Evans,^'^'^^ Huan Lin,'^''^'^ and F. D. A. 
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ABSTRACT 

We have previously identified an excess population of predominantly red galaxies around a 
sample of 31 radio-loud quasars (RLQs) at 1 < z < 2. Here we show that these fields have a 
surface density of extremely red objects (EROs, with R — K > 6) 2.7 times higher than the 
general field. Assuming these EROs are passively evolved galaxies at the quasar redshifts, they 
have characteristic luminosities of only ~ L* . Only one of four RLQ fields has an excess of 
J — K selected EROs with J — K > 2.5; thus, those objects are mostly unrelated to the quasars. 

We also present new multiwavelength data and analyses on the fields of four of these quasars 
at Zq 1.54, obtained to build more detailed pictures of the environments of these quasars and 
the galaxies within them. First, wide-field J and Ks data shows that the galaxy excess around 
Q 0835-1-580 is of Abell richness 2±1 and extends to 140", and that the galaxy excess around 
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Q 1126+101 extends to only 50" civcii though the overall counts in the field are higher than the 
literature average. Second, in three fields wc present the deepest narrow-band redshifted Ha 
observations yet published. We detect five candidate galaxies at the quasar redshifts, a surface 
density 2.5 times higher than in the only existing random-field survey of similar depth. However, 
photometric spectral energy distribution (SED) fitting of one candidate suggests it is an [O III] 
detection background to the quasar. Third, SCUBA submillimeter observations of three fields 
detect two of the quasars and two galaxies with SEDs best fit as highly reddened galaxies at the 
quasar redshifts. Fourth, H-h&nd adaptive optics (AO) imaging is used to estimate redshifts 
for two moderately red bulge-dominated galaxies in the Q 0835-1-580 field using the Kormendy 
relation between central surface brightness and half-light radius. Both have structural redshifts 
consistent with early-typo galaxies foreground to the quasar at z < 0.2 or 1 < z < 1.35. 
Photometric redshifts do not confirm these structural redshifts, however, possibly because our 
optical photometry for those objects is corrupted by scattered light from the nearby bright AO 
guidestar. Finally, quantitative SED fits are presented for numerous galaxies of interest in two 
fields and are used to constrain their photometric redshifts Zph- Most galaxies in the spatially 
compact group around Q 0835-1-580 are consistent with being at the quasar redshift Zq. One of 
these is a candidate very old galaxy without ongoing star formation, while the others appear 
to have ongoing or recent star formation. Many very and extremely red objects across both 
fields have Zph ~ Zq, and significant dust is required to fit most of them, including about half 
of the objects whose fits also require relatively old stellar populations. Large reddenings of 
E{B - y) ~ 0.6 ± 0.3 are also required to fit four J - K selected EROs in the Q 1126+101 
field, though all but one of them have best- fit redshifts Zph > Zq. These objects may represent 
a population of dusty high-redshift galaxies underrepresented in optically selected samples. 

Taken together, these observations reinforce the claim that radio-loud qiiasars at Zq > 1 can 
be found in galaxy overdensities. Ongoing star formation with moderate amounts of dust seems 
to be common among all but the very reddest galaxies in these overdensities. 

Subject headings: infrared: galaxies — galaxies: clusters: general — galaxies: General — 
quasars: individual (Q 0835+580, Q 1126+101, Q 2149+212, Q 2345+061) 



1. Introduction 

It is of considerable current interest to eflaciently identify clusters and other concentrations of galaxies at 
z > 1 in order to study the evolution of both galaxies and large scale structure. Radio-loud quasars (RLQs) 
are among the obvious signposts around which to search for clusters at z > 1 . This is because RLQs almost 
exclusively reside in giant elliptical galaxies (a strongly clustered population) and some RLQs have been 
spectroscopically confirmed to inhabit clusters at ^; < 0.7 (Yee & Ellingson 1993). In Hall, Green & Cohen 
(1998), hereafter HGC98, and Hall & Green (1998), hereafter HG98, wc presented optical (r < 25.5) and 
near-IR [K < 21) imaging of the fields of 31 RLQs at 1 < z < 2 which revealed an excess of predominantly 
red galaxies at K > 19. The excess has two spatial components — one within 40" of the quasars and one 
extending out to at least 2' — consistent with the quasars often residing in galaxy groups or poor clusters 
which are themselves found within richer filaments or sheets of galaxies. The images of HGC98 were at most 
3'x3' in size, and so could only place a lower limit on the size of the large-scale galaxy excess (relative to 
blank-field galaxy counts at these magnitudes) around the quasars. The magnitude and color distributions 
of the excess galaxies are consistent with a population of predominantly early-type galaxies at the quasar 
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redshifts. However, a possible exception is some or all of the J-band dropouts {J — K > 2.5) found in five 
fields for which we have J imaging. Such very red J — K colors require that the galaxies are reddened by 
dust, are at z > 2.5 so that the 4000 A break lies beyond the J band, or both. 

Spectroscopic confirmation of overdensities at the quasar redshifts in these RLQ fields is still lacking, as 
are spectroscopic redshifts for the J-band dropouts, due to the difficulty of obtaining spectroscopic redshifts 
for galaxies at 1.5 < z < 2.5, and extremely red galaxies at all ^; > 1. To provide a more comprehensive 
picture of these systems in the absence of spectroscopy, in this paper we present further analyses of the excess 
galaxy population as well as various new observations of four of these RLQ fields. The results reported here 
supersede earlier results reported as work in progress in Hall et al. (1999) and Hall, Sawicki & Lin (2000). 
In §2 we present new wide- field near-IR imaging and study the extent of the galaxy overdensities. In §3 we 
discuss the surface density of extremely red objects (EROs) in these RLQ fields compared to field surveys. 
In §4 we present narrow-band imaging sensitive to Ha emission from galaxies at the quasar redshifts. In 
§5 we present sub-millimeter mapping sensitive to dust in luminous star-forming galaxies. In §6 we present 
adaptive-optics imaging of two galaxies. In §7 we present multicolor SED fits used to derive photometric 
redshifts, reddenings, and age estimates for selected galaxies. We summarize our results in §8. We adopt 
a A = cosmology with Ho=75 km s~^ Mpc"-^ and go=0.1 (for a projected scale of 7.4 /175 kpc/arcsec at 
z = 1.54). except as noted for comparison with the literature. 



2. Wide-Field Near-IR Imaging 

We obtained J and Kg data over wider fields around two RLQs from HGC98 to help verify the reality 
of the large-scale galaxy excess and to search for additional J-band dropouts. Information about these 
RLQs, the two other RLQs studied in this paper, and the observations obtained in the field of each is 
given in Table 1. The near-IR imaging data were reduced using PHIIRS.^'^ Standard infrared observing and 
reduction procedures were used, including flattening using either domeflats or "running skyflats" created 
from the science images themselves. Object detection and photometry using FOCAS (Valdes 1982) followed 
the procedure of HGC98. 



2.1. Q 0835+580 Field Observations and Reductions 

We obtained J and Ks data over a ~9' diameter unvignetted field of view at platcscale O'.'5/pixel at the 
Steward Observatory 90" telescope with PISCES (McCarthy et al. 1998, 2000), which incorporates a 1024^ 
HgCdTe detector. PISCES has significant optical distortions in both J and Kg which must be removed 
before coadding. Observations of M 67 (Girard et al. 1989) were used to calculate distortion solutions for 
the April 1999 run, yielding images oriented to true North and East (as defined by the USNOA-2.0 catalog) 
to within <0.12° accuracy. Calibration used ARNICA (Hunt et al. 1998) and NICMOS (Persson et al. 
1998) standards and assumed extinction coefficients of 0.08 in J and 0.07 in Kg. Both object detection and 



"PHIIRS is Pat Hall's Infrared Imaging Reduction Software (Hall, Green & Cohen 1998), a package of IRAF 

routines available from the first author or the IRAF contributed software website. The Image Reduction and Analysis 
Facility (IRAF) is distributed by the National Optical Astronomy Observatories, which is operated by AURA, Inc., 
under contract to the National Science Foundation. 
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aperture definition were done on the summed J + Kg image (total area 52.561 arcmin^), and FOCAS total 
magnitudes were produced using those same apertures on the J and Kg images. 



2.2. Q 1126+101 Field Observations and Reductions 

We imaged Q 1126+101 with the near-infrared imager /spectrograph TIFKAM (a.k.a. ONIS at Kitt 
Peak; Pogge et al. 1998) at the 2.4m Hiltner telescope of the MDM Observatory. TIFKAM has a 512 x 1024 
InSb array with a platcscalc of O'.'3/pixcl at the 2.4m, corresponding to a 2.5' x 5' field of view. We observed 
a field centered 80" W of the quasar (which includes the quasar near the eastern edge of the array) in J 
and Ks in nonphotometric conditions of ~1'.'5 seeing. The coaddcd images were trimmed to the region of 
maximum exposure time (12.051 arcmin^). The final J and Kg images completely overlap the r— band data 
of HGC98 and partially overlap the J and Kg data of HGC98. We rotated and magnified our new near- 
infrared observations to match those images and then catalogued all the objects in the sum of the r image 
and our new J and Kg images. Finally, we calibrated our J and Kg images using the objects in common 
with HGC98. 



2.3. Wide-Field Near-IR Imaging Results 

The K < 20.5 galaxy surface density is 17.5±0.6 arcmin'^ in the Q 0835+580 field and 20.7+1.3 
arcmin"^ in the Q 1126+101 field. These values are respectively 2.4o- and 3.5a above the literature com- 
pilation value of 13.7+1.5 (including RMS field-to-field scatter) from HG98. Both values agree within 1.3a 
with those measured in HGC98 using smaller (~3'x3') images of these fields. The J < 21.5 galaxy surface 
density is 11.8+0.5 arcmin'^ in the Q 0835+580 field and 14.2+1.1 arcmin^^ ^hc Q 1126+101 field. These 
values are slightly higher than the J-band field galaxy counts of Teplitz, McLean & Malkan (1999) and > 3a 
higher than those of Saracco et al. (1999). Thus these new wide-field observations confirm that the surface 
density of faint J- and /C-selected galaxies is higher in these RLQ fields than in random fields. 

The radial distribution of all K < 20.5 galaxies relative to Q 0835+580 (Figure la) shows a clear excess 
at <20". There also appears to be an overdensity extending out to ^140" (^ "^h^^ Mpc). From the surface 
density at >140", we would expect 245 galaxies within 140" of the quasar. We observe 285, which is a 2.35a 
excess. The radial distribution of all K < 20.5 galaxies relative to Q 1126+101 (Figure lb) shows a constant 
surface density within the uncertainties to 180", apart from a weak excess at < 50". Given this relatively 
flat radial distribution, part of the excess population of faint J- and X-sclcctcd galaxies in the Q 1126+101 
field is probably unrelated to the quasar. This is consistent with the J — K and r — K color distributions, 
which show an excess of blue galaxies as well as of the red galaxies we expect to find at the quasar redshifts. 
Another possibility is a structure at the quasar redshift of size > Qh^^ Mpc (cf. Sanchez & Gonzalez-Serrano 
1999), which is less probable but cannot be ruled out without spectroscopy or data on even wider fields. 

As a measure of the richness of the structure around Q 0835+580, we computed A^o.s (Hill & Lilly 
1991), the excess number of galaxies within 0.5 Mpc of the quasar and no more than 3 magnitudes fainter 
than a brightest cluster galaxy at the quasar redshift. A^o.5 was computed as in HG98 except that galaxies 
>140" from the quasar were used as the background, rather than the average literature counts. We find 
A^o. 5=27+11 (corresponding to Abell richness 2+1) compared to the value of 16+13 found in HG98. This 
newer, more robust measurement illustrates the value of wide-field data in determining the background 
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counts locally rather than taking them from the literature or from small radii which may still be inside the 
large-scale structure being studied. In the case of Q 1126+101, even though the overall counts in the field are 
higher than the literature average, wide-field data suggests that the quasar is only embedded in a small-scale 
overdensity. 



3. Surface Density of EROs in 1< z < 2 RLQ Fields 

Since the appearance of infrared array detectors, numerous discoveries of extremely red objects (EROs) 
have been noted in the literature (for recent summaries see Dcy et al. 1999; Daddi et al. 2000). EROs remain 
a poorly understood population of galaxies about which several different interpretations are plausible. Their 
observed colors oi R — K > & are so red that they seem to be explainable only by old stellar populations 
at z > 1 or by heavily dust-reddened galaxies or AGN preferentially located at z > 1. This is confirmed by 
the handful of spectroscopic redshifts for EROs (Dey et al. 1999; Soifer et al. 1999; Liu et al. 2000) and very 
red objects (VROs), as defined by (Cohen et al. 1999, ; cf. §7.3). Previous narrow-field optical and near-IR 
surveys have hinted that EROs are more common along lines of sight to distant AGN (for recent summaries 
sec Small et al. 1999; Cimatti et al. 2000). Our 1 < z < 2 RLQ fields can be used to test this hypothesis 
more carefully, using recent estimates of the surface density of field EROs from wide-field surveys. 



3.1. r-K selected EROs 

Thompson et al. (1999b) find 29 K' < 19 EROs with R - K' > Q m four fields of ~150 arcmin^ 
each, for a surface density of 0.048^g Q|!)g arcmin"^ (we use the methods of Gehrels (1986) to find all la 
Poissonian uncertainties where small number statistics apply). The details for one field have been published 
in Thompson et al. (1999a; hereafter T99). The T99 K' observations were calibrated using UKIRT standards 
(Gasali & Hawarden 1992), as were the Kg observations of HG98. We refer to both magnitudes simply as 
K hereafter. We also adopt R — K > Q ds the definition of a ERO, or r — > 6.322 for Gunn r imaging 
such as ours (HG98, Appendix A), and X = 19 as the division between bright and faint EROs. However, 
we conservatively adopt a 3(7 detection limit in the optical instead of 2a as in T99. This yields a difference 
of ~0™44 in the lower limit on the r — K color of undetected objects. 

We select EROs from all 22 RLQ fields of HGC98 with 3a limits of r > 25.322 and 5a limits of 
if > 19 and from the new larger area aroimd Q 1126-1-101 imaged with TIFKAM (§2.2). The total area is 
160.95 arcmin^, almost identical to that of T99, and the depth of our imaging is also similar. We find 24 
EROs, including two presumed stars unresolved in archival HST WFPC2 snapshots and one presumed star 
identified via its blue J — K color. Four of the remaining 21 EROs are compact enough that they could be 
stars given the seeing, but the rest are extended. T99 found two of their eight bright EROs to be unresolved 
in I'.'l seeing, and both were subsequently confirmed to be stars. There is no reason to expect more stellar 
contamination in our ERO sample than in the T99 sample, since both are from high Galactic latitude fields 
of almost exactly the same total areas. Therefore we assume that only three of our 24 EROs are stars. The 
remaining 21 r — K selected extragalactic EROs yield a surface density of 0.130tQ Q29 arcmin^^ to K = 19, 
which is 2.7 times higher than the field ERO counts of Thompson et al. (1999b) and a 2.64cr excess over 
the number expected in this area based on those counts. This is consistent with the excesses found by the 
smaller surveys of Chapman, McCarthy & Persson (2000) at if < 20 around 4 radio galaxies with z ^ 1 and 
Cimatti et al. (2000) at < 19 in 40 arcmin^ around 14 radio-loud AGN with z > 1.5. 
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3.2. J -K selected EROs 

Drory et al. (1999) find a surface density of 0.098±0.009 arcmin-^ for K < 19.5 EROs with J-K>2.^ 
from 124 objects over 1260 arcmin^ from the random-field MUNICS survey (Mendes de Oliveira et al. 1998). 

Stellar contamination of these J — K selected EROs should be negligible. A literature search revealed no 
predicted colors oi J — K > 2.2 for unreddened stars of any type from brown dwarf to supergiant and very 
few observations of objects with J — K > 2.5 which were not galaxies: two objects presumed to be stars 
since they were unresolved on HST WFPC2 images (De Propris et al. 1999) and a small number of known 
AGN (Cutri et al. 2000; Barkhouse & Hall 2001). Thus objects with J - K > 2.5 are almost certainly not 
stars, are unlikely to be AGN, and are probably galaxies. 

In four RLQ fields from HGC98 with deep J data (Q 0835+580, Q 0952+179, Q 1126+101, and 
Q 1258+404), we find six such J — K selected EROs. We also find three J — K selected EROs in the new 
area around Q 0835+580 imaged with PISCES, but none in the new area around Q 1126+101 imaged with 
TIFKAM. Five of the EROs are from the Q 1126+101 field, which as discussed in HG98 has a considerable 
excess of these J-band dropouts. Excluding the Q 1126+101 field because of its excess number of these 
objects, we find four J — K selected EROs in 70.373 arcmin^ for a surface density of 0.057lo!o24 arcmin"^, 
within la of the field measurement of Drory et al. (1999). In the Q 1126+101 field we find 3 EROs in 
15.549 arcmin^ for a surface density of 0.19^Q jg arcmin"^, 3.3 times higher than in random fields but only 
significant at 1.54(t due to small number statistics. 

We also consider the surface density of J-band dropouts in the Q 1126+101 field to K = 20, since 
all our data reaches deep enough to select EROs with J — K > 2.5 to at least this magnitude. In the 
new TIFKAM Q 1126+101 data we find 6 such EROs (0.50to:20 arcmin^^). In the HGC98 Q 1126+101 
data we find 8 such EROs (0.93^032 arcmin^^). In the other three fields with deep J data from HGC98 
plus the PISCES Q 0835+580 data we find 11 such EROs (0.16i[J;JJ^ arcmin'^). The enhancement in the 
Q 1126+101 field is still a factor of three or more, but remains only 2(j significance due to small number 
statistics. The same holds true down to K = 20.5 in the HGC98 data alone. Deeper or wider-field data 
are unlikely to dramatically increase the statistical significance of the overdensity of J-band dropouts, but 
their clustered distribution within the Q 1126+101 field (Figure 24 of IIG98) and the overdensity of galaxies 
with less extreme colors suggests that this is a physical concentration and not just a random fluctuation. 
However, some recent work suggests that such random ficld-to-field fluctuations can be quite substantial 
even to K = 2Q and beyond (cf. Daddi et al. 2000; Eisenhardt et al. 2000; Scodeggio & Silva 2000). 



3.3. Extremely Red Objects: Discussion 

We find bright r — K selected EROs with < 19 to be 2.7 times more common in the fields of 22 RLQs 
with 1 < < 2 than in the general field. We find J — K selected EROs with K < 19.5 to be from 1 to 3.3 
times as common in the fields of four RLQs with 1.5 < Zq < 1.75, depending on whether the Q 1126+101 



^''Two of these EROs in the Q 1126+101 field have inconsistent measurements (at > Scr) between the old HGC98 
and new TIFKAM datasets. Examination of the images suggests that their J fluxes were probably underestimated 
in HGC98. The HGC98 KPNO 4m TRIM data required 'destriping' to remove the large scale scattered light pattern 
present in some fields (see §3.2.7 of HGC98). The new J data shows that this procedure worked quite well overall, 
but in a few cases the masking of objects was insufficient to prevent removal of some real flux. 
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field is included or not. These quantitative measurements show that the excess ERO surface density around 
distant AGN is considerably smaller than previous estimates of a factor of 10 to 100 excess (Dey, Spinrad & 
Dickinson 1995). 

One outstanding question regarding EROs is their intrinsic luminosities. The presence of the quasars 
in these fields enables us to estimate the luminosities of the r — K selected EROs by assuming that the 
EROs are at the quasar rcdshifts. This is a reasonable assumption since there is no evidence for strong or 
weak gravitational Icnsing in these fields and since the excess galaxy population is concentrated around the 
quasars but not correlated with the presence of intervening Mgii absorption systems (see §2 and §4.3 of 
HG98 for details and discussion, and §7.3 for further evidence of association). 

Assuming the EROs in each field are at the redshift of the quasar in that field (<Zg>=1.61) and applying 
elliptical fc-corrections only from Fioc & Rocca-Volmerange (1997), we find M^^'^=-25.4±0.5 for the r — K 
selected EROs in our standard cosmology. Applying evolutionary corrections (e-corrections) from Poggianti 
(1997) as well, we find M^-^'^=-24.1±0.4. This is consistent with the Gardner et al. (1997) measurement 
of M^=-24.2±0.2 for field galaxies at <2;>=0.14 (see HG98 §5.1) and with the De Propris et al. (1998) 
measurement of M^=-23.9 in the Coma cluster. 

Thus, these EROs have absolute magnitudes consistent with those expected for passively evolving 
luminous elliptical galaxies at the quasar rcdshifts, and need not be an extraordinarily luminous galaxy 
population. The higher density of EROs in RLQ fields can be easily understood if the RLQs are located in 
overdense regions. Our measured ERO luminosities are consistent with the estimated luminosity evolution 
in for excess galaxies to > 20 in these fields (§5 of HG98; see also Kajisawa et al. 2000). The faint 
end of our ERO distribution is set by our K <19 selection criterion, however, so we can say nothing further 
about the ERO luminosity function. 

Using k- and e- corrections for bluer galaxies of later spectral types yields higher estimated luminosities 
for the EROs, but the red colors of these EROs strongly suggest that such corrections are not appropriate. On 
the other hand, if dust reddening is important for these EROs at 2; > 1 but not for their descendants at 2; = 0, 
strong k- and e- corrections would be required. If dust were important we might expect considerable overlap 
between r-K and J — K selected EROs. Thompson et al. (1999b) find that only A oi 29 r — K selected EROs 
are J — K selected EROs, and of our four objects with r — K> 6.322 and J data, only one has J — K > 2.5. 
Our SED fitting (§7.3) yields a best-fit E{B - l/)=0.75±0.50 for this object, E{B - y)=0.45±0.20 for one of 
the J — K < 2.5 objects, and E(B — V)=0 for the remaining two J — K < 2.5 objects. If the same fractions 
hold among the r — K selected EROs of Thompson et al. (1999b), up to ~12 of 29 (~40%) could have strong 
dust reddening. Multicolor data on larger r — K selected ERO samples reaching fainter magnitudes will be 
needed to determine how frequent and strong dust reddening really is among EROs, but for now we make 
no correction for dust reddening. 



4. Narrow-band Ha Imaging 

Ha is a good tracer of the instantaneous star formation rate (SFR) in galaxies (Kennicutt 1998). 
Narrow-band surveys until recently yielded only limits on the space density of Ha emitters at 2: > 1, even in 
the fields of quasars and radio galaxies (Pahre & Djorgovski 1995). However, in the past few years detections 
have been made of objects located in random fields and objects associated with damped Lya absorbers, with 
strong intervening metal-line absorbers, or with quasars or radio galaxies (for a recent overview see Teplitz, 
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McLean & Malkan 1999). Given the excess galaxy population scon in our RLQ fields, they make promising 
targets for narrow-band observations to search for Ha emission at the quasar redshifts. In particular, we 
would hope to detect galaxies whose SEDs suggest they are dust-reddened and thus possibly actively star- 
forming and we would expect not to detect galaxies whose SEDs suggest they are old and dust-free. 

The fields of Q 0835-^580 and Q 2345-^061 were imaged in H at O'.'300/pixel and in a narrow band 

within H at O'.'148/pixel using NSFCAM (Leggett & Deanult 1996) at the NASA Infra-Red Telescope Facility 
(IRTF) . Relevant details of the observations are given in Table 1 . Narrow-band imaging utilized a circularly 
variable filter (CVF) with resolution R = 90 tuned to the wavelength of Ha at the quasar reshifts. Photo- 
metric calibration was done using UKIRT standards (Casali & Hawarden 1992). Domeflats were used in H 
but not for CVF imaging since the fringing across the array in CVF mode was different between domeflats 
and the sky. Photometry was performed in matched 4.2" (14 pixel) diameter apertures on the H images 
and on Ha images resampled to the i/-band pixel scale. Additional H- and narrow-band observations of the 
fields of Q 2345-1-061 and Q 2149-1-212 were obtained by C. Pritchet & F. Hartwick (Table 1). These data 
were taken with CFHT and REDEYE (O'.'5/pixel) and a custom-made narrow-band filter 62 A wide centered 
at 1.6662 /xm. The data were calibrated using bright objects in common with the IRTF Q 2345-1-061 data, 
and magnitudes were measured in 4'.'0 (8 pixel) diameter apertures. To quantify the significance of the excess 
narrow-band flux we follow Bunker et al. (1995) and use S, the number of standard deviations between the 
observed broad-band counts and the number expected based on the narrow-band counts. 

There is one > 3a detection in each field observed with IRTF (see Table 2). Figure 2a shows the 
color-magnitude diagram for the Q 0835-1-580 field. The Ha emitter — hereafter Q 0835-1-580 (Hal) — is 
an unremarkable faint blue galaxy with SFRHa=14.7±2.5 h^^MQyr~^ using the Kennicutt (1983) relation 
between SFRho and L(Ha), which has an additional intrinsic dispersion of ±50%, and assuming the detection 
is indeed Ha at the quasar rcdshift. With this assumption, we can also estimate SFRi?t/v following Steidel 
et al. (1996) using i7-band data which has been obtained for this field (see §7), since U samples rest-frame 
1500 A at the quasar redshift almost exactly. Q 0835-1-580 (Hal) has ?7=23.052±0.075, which corresponds to 
SFRir(7v=5.3±0.3 h^^MQyr~^ for a Salpeter IMF with an upper mass cutoff of 80 Mq. The agreement with 
SFRna is good given the various imccrtainties (e.g. no correction for dust has been made to SFR^-f/y). Note 
that none of the dozen or so galaxies with r — K > 5 immediately surrounding Q 0835-1-580 were detected 
in Ha. If they are at the quasar redshift, they must have star formation rates of < 1.8h^^MQyi~^ (^c), 
consistent with them being red due to age and/or metallicity rather than being extremely dust-reddened 
starburst galaxies (see §7). 

Figure 2b shows the IRTF data color-magnitude diagram of the Q 2345+061 field, and Figure 3a the 
same diagram for the CFHT data. Figure 4a shows the Ha image with the Ha emitters marked. The 
candidate Ha emitter seen with IRTF — Q 2345-1-061 (Hal) — is not confirmed with CFHT, but the CFHT 
narrow-band filter is much narrower than the IRTF CVF and the line could lie outside the wavelength range 
of the CFHT filter. There are also two objects visible in the narrow-band CFHT image but not on the 
if-band CFHT image. They both have Ha excesses of 3a significance with H — Ha > 1.7 (Figure 3a and 
Table 2). Q 2345-^061 (Ha2) is detected at if=20.3±0.1 and H - i?a=0.4±0.3 in the IRTF data while 
Q 2345+061 (Ha3) is not detected in the IRTF Ha image or in the deep data of HGC98 to Ks>21.7. 
It may be detected at ir~21.5±0.5 in the IRTF data and at r~25.1±0.5 in the r data of HGC98. In any 
case, given the photometric uncertainties, the CFHT and IRTF Ha measurements are consistent for both 
Q 2345+061 (Ha2) and Q 2345+061 (Ha3). 

Figure 3b shows the color-magnitude diagram for the CFHT data on the Q 2149+212 field. There is 
one candidate Ha emitter, with H ~ 20.2 (Figure 4b and Table 2). It is a fairly red galaxy with r — K = 5 
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and H — K ^ 1. This object is confirnicd as an Ha emitter in additional data taken the same night but not 
used due to cahbration problems, whereas the object at if ~ 19 and H — Ha 0.3 is not confirmed. 

Our Ha imaging reaches flux limits of 1 to 2 x 10^^'' erg cm^^ s^^ (Sc), deeper than any previously 
published survey due to the very long IRTF CVF exposures and the very low background in the CFHT 
narrow-band images. Five detections in 10.156 arcmin^ (counting both observations of the Q 2345+061 field 
since different filters were used) gives a surface density of 0.5to'2 arcmin"^. For comparison, no Ha sources 
were found in two random fields surveyed to 0.8 x 10~^^ erg cm~^ (3(t) with the same CFHT setup as 
in this work, for a la upper limit of 0.2 sources per arcmin^ (Pritchet et al., in preparation). In the deepest 
published survey, McCarthy et al. (1999) find 0.57 sources per arcmin^ to 4.1 10~^^ erg cm~^ s~^ {3a) from 
slitless spectroscopy of random fields with NICMOS. We find only 0.2 sources per arcmin^ to the same limit, 
but they sample a rcdshift range of Az = 1.2 compared to our average Az = 0.013. Our surface density is 
approximately 32 times higher, although that number is biased high because we targeted a redshift where 
we had reason to believe a galaxy excess existed. In surveys at the redshifts of quasars and radio galaxies, 
Pahrc & Djorgovski (1995), van dcr Wcrf (1997) and Tcplitz, McLean & Malkan (1999) foimd 4 sources in 
242 arcmin^ to limits of 9 to 19 10~^^ erg cm~^ s~^ (^c). We find at most one source in 10.156 arcmin^ to 
these depths, which is consistent with the previous surveys due to small number statistics even though it is 
formally up to a factor of six excess. 

In summary, our fields represent a factor of 2.5 excess compared to the only existing random- field survey 
to equal depth, a factor of <32 excess compared to the deepest published random-field survey, and a factor 
of <6 excess compared to previous AGN-field surveys. The fields studied in this paper are among the best 
candidates for 1 < 2; < 2 RLQs in clusters from a survey of 31 (Hall, Green & Cohen 1998), so this excess is 
not unexpected. These Ha detections constitute more evidence that the galaxy overdensities in these RLQ 
fields are real, and that at least some of the excess galaxies are at the quasar redshifts. Although the number 
of member galaxies in the candidate host clusters of Q 2149+212 and Q 2345+061 is unknown, the deep 
CFHT Ha images show that there are only three such galaxies with star formation rates of > l.bh^^MQ yr~^ 
within fields 0.95h^^ Mpc wide centered on these quasars. This is a lower limit which neglects extinction 
and the velocity dispersion of the clusters (e.g. Q 2345+061 (Hal) is not detected in the narrow CFHT 
filter), but it is still a stringent limit which illustrates the potential of deep wide- field narrow-band data in 
studying star formation rates in high redshift clusters. 



5. Sub- millimeter Mapping 

In the last few years it has become clear that much of the star formation activity in the universe is 
obscured by dust and can be easily detected only at far-IR and sub-mm wavelengths (e.g. Dwek et al. 1998; 
Small et al. 1998). The presence of a number of galaxies with SEDs strongly indicative of substantial dust 
reddening in our RLQ fields suggested that they might be detectable sub-mm sources. Thus we obtained 
continuum observations at 450 and 850/Ltm with SCUBA, the Sub-mm Common-User Bolometer Array (Hol- 
land et al. 1998), on the James Clerk Maxwell Telescope. Details of the observations are given in Table 1. 
All observations used standard 64-point jiggle maps to fully sample both arrays. 

Data on Q 2345+061 were obtained (and reduced) as a service observing program. Azimuthal chopping 
and nodding was used. No sources are seen at 450/im, but at 850/im Q 2345+061 is detected at 2.8a (Table 3). 
No other potential sources (e.g. the Ha emitters) are seen, but we are only sensitive to hyperluminous IR 
galaxies [LpiR > lO^^h'^Q Lq). The 2a limit at 850/im corresponds to ~ 1.3 x lO^^h^^ Lq for the quasar 
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redshift and a bolometric correction calculated using a dust emissivity spectral index of 1.0 and the Arp 220 
dust temperature T = 47 K, following Barger et al. (1999). 

Observations of Q 0835+580 and Q 1126+101 were made using a fixed chop throw of 40" and 12'.'36 
offsets N or S every other measurement. Conditions ranged from fair to excellent during the run. The data 
were reduced and calibrated in the standard manner using SURF (Jenness & Lightfoot 1998). Figures 5 and 
6 show 850/im contours superimposed on K images of the Q 0835+580 and Q 1126+101 fields, respectively. 
The 450/im maps arc not shown as no objects were detected at that wavelength in cither field. To help 
reject spurious sources, we cross-correlated the maps with the beam as given by a calibration observation 
of an unresolved source on the same night. Each pixel's value in the cross-correlation map measures how 
closely the sub-mm map centered on that pixel resembles the beam, providing a measure of the correlation 
coefficient (CC) for each tentative source. Real high-redshift objects should resemble the beam quite closely, 
since objects more extended than a point source cannot be at high redshift given the large JCMT beam and 
objects more compact than a point source might be just noise. 

There are two > 3.5cr sources in the Q 0835+580 field which coincide with peaks in the cross-correlation 
map, but the bright source (SMI; 10.1 mJy) has only CC=0.37 and the faint source (SM2; 7.5 mJy) only 
CC=0.29. There are nearly equal numbers of negative and positive peaks at these fiux and CC values, so 
we regard both these sources as tentative sub-mm detections. However, they both have possible optical/IR 
IDs: aK = 19.8 moderately blue galaxy (Q 0835+580 (117), r - K=Z.Q) 3'.'2 from the position of SMI, and 
a.K ^ 17.7 red galaxy (Q 0835+580 (458), r - i^=5.4) 2'.'1 from the position of SM2. The ID numbers of 
these objects refer to our optical/NIR catalogs of these fields, available on request from the first author. We 
include information on these objects in Table 3 (and Table 4), but we emphasize that we do not consider them 
unambiguous sub-mm detections. Both objects are faint enough in the optical and near-infrared that they 
could be at the quasar redshift; Q 0835+580 (SM2) would be 0™4 brighter than the magnitude of a brightest 
cluster galaxy (as estimated in HG98) if it was at the quasar redshift, but this is within the observed scatter 
in such objects' magnitudes at ~ 1 (Thompson, Aftreth & Soifer 2000). We discuss these objects further 
in §7.1.4. Note that Q 0835+580 (Hal) is not detected, although the 3a limit of SFRpiR < 2090 ± 550 
K^^Mq yr~^ (calculated following Dey et al. (1999) using a Salpcter IMF and assuming it is at the quasar 
redshift) is two orders of magnitude larger than SFRf/y and SFRff^ for this object (§4). 

In the Q 1126+101 field, the quasar and several other sources are present with > 3cr and CC>0.33. 
There are twice as many fluctuations above sky as below sky at these flux and CC values in this field; 
thus, some of the detections are probably real but some are probably spurious. Overall the most reliable 

detections are the quasar and SMI. Despite their low S/N, they were clearly detected on both nights, have 
high correlation coefficients, and have FWHM~14" as expected for the JCMT beam at 850/xm. The possible 
counterpart I'.'l to the N of SMI, object Q 1126+101 (384), is discussed further in §7.2.1. The two other 
sources listed in Table 3 are less certain, despite their higher S/N on the raw map. They have FWHM^9", 
suggesting they might be noise spikes rather than actual objects. Neither Q 1126+101 (SM2) nor (SM3) 
have probable IDs in the optical/IR data. 



6. CFHT Adaptive-Optics Imaging 

Little information is available on the detailed morphologies of EROs, although our knowledge is growing 
rapidly (Dey et al. 1999; Thompson et al. 1999a; Benftez et al. 1999; Moriondo, Cimatti & Daddi 2000; Sti- 
avelli & Treu 2000). To improve this situation, il-band adaptive-optics data were obtained with the CFHT 
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Adaptive Optics Bonnette (known as Pue'o) and the KIR camera (1024^ HgCdTe detector, O'.'0348/pixel). 
The y=10.3 star BD+58 1138 was used to guide on a field -35" from the star and ~90" from Q 0835+580. 
Standard infrared observing and reduction procedures were used, including domeflats for flattening. Photom- 
etry was done using UKIRT standards and assuming a typical Mauna Kea H-h&nd extinction coefBcient of 
0.051±0.012 (Krisciunas et al. 1987). Natural seeing conditions during the observations were poor, resulting 
in FWHM after correction of ~ 0'.'42, as measured from the most compact object on the image. 

The KIR field contains two objects bright enough for image profile analysis. The fainter, redder object 
is designated Q 0835+580 (106) since it is #106 in our catalog of this field and has ^8=17. 7, r — K=^ 
and J - K=2. The other object, Q 0835+580 (112), has Ks=n.^ and r - K=2A. Radial profile analysis 
shows that both galaxies arc much better fit by an r^/^-law profile than by an exponential disk profile, with 
re=0'.'97±0'.'14 for Q 0835+580 (106) and re=0'.'49±0'.'07 for Q 0835+580 (112). 

The resolution was not high enough for further morphological analysis, but if we assume the objects 
are early-type galaxies as indicated by their radial profiles, we can use the Kormendy relation between 
effective radius and average effective surface brightness (Kormendy & Djorgovski 1989) to determine what 
redshifts are consistent with the objects' observed sizes and surface brightnesses (Eisenhardt et al. 1996). 
To perform this analysis we use the elliptical colors and evolutionary and fc-corrections of Poggianti (1997) 
and adopt their iJo=50, (70=0.225 cosmology for consistency. For each redshift at which Poggianti (1997) 
calculate e- and fc-corrections, we use the CFHT data to compute the rest-frame B-band surface brightness 
within (accounting for (l+z)^ surface brightness dimming) and the physical size corresponding to the 
observed r^- The tracks for both objects are shown in Figure 7 along with data for local ellipticals from 
Sandage & Perelmuter (1990). Both objects intersect the blue surface brightness-effective radius relation for 
ellipticals both at low and high redshift. Q 0835+580 (112) could be at 0.04 <z< 0.24 or 1< ^ < 1.4, and 
Q 0835+580 (106) could be at 0.02 < z < 0.10 or 1 < < 1.3. We compare these estimates with those from 
photometric redshifts in §7.1.3. 



7. Photometric Redshifts and Other Constraints from SED Fits 

Estimating photometric redshifts Zph and spectral types from multicolor photometry can provide insights 
into the nature of faint, high- 2; objects which are otherwise too faint or too numerous for efficient spectroscopy 

(for a recent review sec Yce 1998). In our RLQ fields we wish to quantitatively confirm or refute the 
existence of excess galaxies at or beyond the quasar redshifts and to investigate the SEDs — and hence the 
evolutionary state — of such galaxies, particularly EROs. A longer-term goal is to provide an efScient way 
to reject contaminants foreground to the quasars in future spectroscopic followup studies. Currently, we 
have sufficient multicolor imaging data for such fitting in only the Q 0835+580 and Q 1126+101 fields. 

An established technique for measuring photometric redshifts is to compare the photometry of objects 
of interest with that of an empirical template set. This technique works well for galaxies with normal SEDs, 
where local template spectra (e.g. the set of Coleman, Wu, & Weedman 1980) provide good results even to 
high redshifts (e.g. Sawicki, Lin & Yee 1997). However, because EROs are atypical objects, and potentially 
very dusty, it is inadvisible to fit their photometry with templates that are representative only of normal 
galaxies. Instead, it is important to allow for the possibility of heavy dust obscuration using synthetic 
templates reddened with a range of dust values. 

We thus constructed a grid of models spanning a range in redshift, age, and reddening. The basis of 
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this grid is the Bruzual & Chariot (1996) version of the spectral synthesis models of Bruzual & Chariot 
(1993), spanning 221 age steps < i < 2 x 10^° yr. We considered instantaneous burst and continuous 
star formation models, which bracket the extremes of monotonically decreasing star formation histories. We 
assumed a metallicity of 0.2Zq, though solar metallicity models were also considered for a subset of objects. 
Our synthetic spectra were attenuated using the Calzetti (1997) starburst galaxy extinction recipe, for a 
total of 67 values of E(B — V) spanning < E{B — V)<2. The reddened spectra were then redshifted onto 
a grid of 101 redshifts spanning < z < 5, after which they were convolved with filter transmission curves 
to compute the final grid of galaxy colors as a function of redshift, reddening, and age. 

The observed photometry for each object was compared with the model grids and the best-fitting 

position in Zph, age, E{B — V) space was found by means of a maximum likelihood test. Almost all the input 
data used in these calculations was reported in HGC98. The exception is five hours of U imaging of the 
Q 0835+580 field obtained with the Steward Observatory 90" on UT 1997 December 31 and reduced in the 
same manner as the data of HGC98. FOCAS total magnitudes were measured in all filters using apertures 
from the summed r + J + Kg image and converted to fluxes using the conversions in Appendix A of HGC98 
plus that of Fukugita, Shimasaku & Ichikawa (1995) for the U data. 

Our procedure is thus similar to that used by Sawicki & Yec (1998) in their investigation of the spectral 
energy distributions of Lyman break galaxies. However, unlike in Sawicki & Yee (1998) where the redshifts 
of objects of interest were known a priori, here we must find the best fitting model in the three-dimensional 
space that covers redshift in addition to age and dust attenuation values. 

For each object fit, xi (reduced x^) contour plots with projected 90% confidence error contours in the 
z/E{B ~ V) and a.ge/ E{B — V) planes were constructed for each SFR scenario (Figures 9-14 and Table 4). 
All our quoted uncertainties and plotted error bars are these projected 90% confidence limits. They are thus 
somewhat more conservative that the standard la (68% confidence) gaussian error bars. It should also be 
kept in mind that the 90% confidence regions are typically not just elliptical gaussians, but are elongated 
because of degeneracies between parameters such as age and E{B — V). Plotting simple error bars, while 
convenient, thus tends to exaggerate the apparent uncertainties. We define a successful fit as one where at 
most one of the three variables of interest — Zph, E{B — V) and age — is unconstrained at 90% confidence 
(i.e., that variable has 90% confidence limits equal to its range of permitted values). Note that a successful 
fit can still be a poor fit with a high value, and vice versa. There is a continuum of properties across 
our dividing line between successful and unsuccessful fits, of course, but objects not successfully fit typically 
have large photometric errors which prevent their SEDs from usefully constraining the fitted parameters. 
However, unsuccessful fits can also be caused by star formation histories which deviate strongly from the 
instantaneous burst or continuous histories we consider, such as multiple strong bursts or slowly declining 
star formation rates. 

In the following sections we present our photometric redshift results. In addition to fitting a handful of 
specific objects of interest, we examine a compact clump of predominantly red galaxies around Q 0835+580 
to confirm or refute the visual impression that these galaxies are associated with the quasar. We also examine 
the population of very red galaxies in these fields regardless of their distance from the RLQ, because these 
galaxies may be part of extended structures at the quasar redshift. Furthermore, these red objects may be 
dusty star bursting galaxies or high-redshift galaxies with old stellar populations, both of which are of great 
interest to the study of galaxy evolution. 
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7.1. Q 0835+580 Field 

The field of Q 0835+580 has data in seven filters (UrizJHK) for most objects. As discussed above, we 
do not compute Zph values for all galaxies in the field, but instead focus in the following sections on specific 
galaxies or galaxy populations. 



7.1.1. The Spatially Compact Group Around Q 0835+580 

There are 18 objects within 2T.'5 of Q 0835+580 to K = 21.24 (the 5a limit of our data) which form a 
spatially compact group in the terminology of Stanford et al. (1997). This count excludes the quasar, one 
bright morphologically classified star, and one probable star which is unresolved on an iJST WFPC2 snapshot 
image. As discussed in HG98 (see their Figure 23), the compact spatial distribution of these galaxies and 
the red color of many of them makes it very likely that they are at the redshift of the quasar (zg=1.5358) or 
the intervening Mgll systems (2;~1.4368). 

To test this hypothesis, we have fitted the broadband magnitudes of these 18 objects. They are denoted 
'SCG' in the Notes column of Table 4, which summarizes the fit results. All objects had successful fits, as 

defined in the previous section. Figure 8a shows E{B — V) vs. Zph (at the best fit model age) for the twelve 
objects with K < 2Q (large symbols) and the six with K > 20 (small symbols). Figure 8b shows the best 
fit model age vs. Zp^. Instantaneous bursts (squares) are preferred over constant star formation histories 
(triangles) in nine of twelve cases at if < 20 but only two of six cases at if > 20. Figure 8b shows the best 
fit model age vs. Zp^. The dotted vertical line in both figures shows the quasar redshift. Recall that the 
simple error bars we plot for convenience do exaggerate the uncertainties somewhat. 

Twelve of the 18 objects are consistent with being at the quasar redshift at 90% confidence. Another 
which misses being consistent by A2; = 0.06 is coincident with the southern radio hotspot of the quasar 
(see §7.2 of HG98), so we consider it to be at the quasar redshift as well. Thus there is good evidence 
from photometric redshifts for galaxy overdensities at the quasar redshifts. For these 13 galaxies the typical 
E[B — V) is ~0.3 and the typical fitted age is ~1 Gyr. Both quantities have large uncertainties, and there is 
a large scatter in the ages. Nonetheless, overall the plots are consistent with current or recent star formation 
in the majority of the galaxies at the quasar redshifts, as denoted by a best-fit constant star formation 
history or a best-fit age of < 1 Gyr. 

The notable exception to this is Q 0835+580 (339), which is a bright ERO with r - K >7 best fit by 
an instantaneous burst with Zph = l.est^j^, E{B -V) = 0.00+!^oJ^, and age 19.75;?3^^o Gyr (Figure 9a). 
Four other objects have similar best-fit ages which are older than the universe at the quasar redshift for 
all reasonable cosmologies, but this is the only one for which even the 90% confidence lower limit on the 
age is older than the universe. For reference, note that for the current best-fit flat, low-f2M, non-zero f2A 
cosmology, at z = 1.5 the age of the Universe was {'i.l+_]yl)h^^ Gyr {^m = 0.2 + 0.1, = 1 - ^m)- To 
see if we have overestimated this object's age by underestimating its metallicity, we refitted it using solar 
metallicity models. Both solar metallicity fits have lower than 20% solar fits, though the instantaneous 
burst is again preferred (Figure 9b), with slightly lower values of all three main parameters: Zpt = 1.40j^Q'2o, 
E{B — V)— O.00lg;og, and age 12.5^8 Gyr. The 90% confidence lower limit on the age for even this solar 
metallicity fit is still older than the age of the universe at the quasar redshift, except for flat A-dominated 
cosmologies with VIm < 0.2. Thus, even though higher metallicities may further reduce the best-fit age, 
Q 0835+580 (339) is clearly worthy of spectroscopic study as a candidate very old galaxy at high redshift. 
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Another of the objects within 21'.'5 of the quasar, Q 0835+580 (352), is known to be at z — 0.236 
(Burbidge et al. 1990), which provides a useful check of the accuracy of our photometric rcdshifts. The 
object is confused with a separate red galaxy (Figure 23 of HG98) which contributes negligible flux in Uriz 
but ~0™5 in J and H and ~0™25 in K. Our fit to its corrected photometry yields Zph = 0.35tolo5 ^'^'^ 
E{B — F)=0.16±0.06. This redshift is only 2.3cr above the known redshift, and this E{B — V) value is close 
to the average value of ~0.1 estimated for field galaxies at <2;>=0.3 (Lin et al. 2001). We thus consider our 
fit to be acceptable. 



7.1.2. Q 0835+580 (Hal) 

The object Q 0835+580 (398) is also known as Q 0835+580 (Hal), a candidate Ha emitter at the 
quasar redshift (§4). However, it is best fit as a young Zp^ — 2.15 ± 0.25 galaxy with E{B — y)~0.35±0.20 
under either the single stellar population or constant SFR models, although a similarly reddened object at 
Zph ~ 0.25 is also possible at 90% confidence. We also fit this object using solar metallicity models and found 
that the best-fit redshift did not change significantly, though the best- fit reddening was somewhat lower, 
E{B — y)~0.20±0.20, and only the constant SFR model still permitted a low-redshift {zph ^ 0.2) solution 
at 90% confidence. It is possible that the narrow-band filter detected [Oiii] A5007 dX z = 2.3338 ± 0.0186 
instead of Ha at 2; = 1.5358, but the Zph ~ 0.25 option would be impossible to reconcile with the narrow-band 
detection. Spectroscopy is needed in this case to verify the photometric redshift. 



7.1.3. Q 0835+580 Field Adaptive Optics Targets 

Photometric rcdshifts were computed for Q 0835+580 (106) and (112) for comparison with the structural 
redshifts calculated for these objects from adaptive optics images (§6). The direct SED fits do not agree very 
well with the structural redshifts. Figure 10a shows that Q 0835+580 (106) is best fit as a Zph = 4.50 ± 0.15 
galaxy with E{B — V) 0.23lg;^g. The structural redshift of z ~ 0.8 is excluded at high confidence for 
the models considered. This result supersedes the earlier modelling attempt reported in Hall, Sawicki & Lin 
(2000). Figure 10b shows that Q 0835+580 (112) is not weU fit by any of our simple models {xl ^ 20). 
A composite stellar population may improve the fit, but that may cast some doubt on its morphological 
identification as an early-type galaxy. From comparison of Figures 10-14, as well as similar plots for other 
objects (not shown), it appears that the riz fluxes for these two objects, especially Q 0835+580 (112), are 
unusually bright relative to the U and JHK fluxes, even compared to blue objects. Thus we believe the 
reason for the poor SED fits to these objects is that their riz photometry is affected by scattered light 
from the nearby adaptive optics guidestar, despite^ attempts to account for such light through masking, 
subtracting, and using local backgrounds. Multiwavelength adaptive optics imaging may be necessary to 
obtain accurate photometry for further study of these objects. 



7.1.4. Q 0835+580 Field Possible Sub-mm Counterparts 

Q 0835+580 (117) and (458) were chosen for SED fitting since they are possible counterparts of tentative 
sub-mm detections (§5). Figure 11a shows that Q 0835+580 (117) is best fit as an instantaneous burst with 
E{B — V)= 0^0^00 ^ redshift consistent with that of the quasar {zph = 1.60 ± 0.15). A constant star 
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formation history requires E{B - V)= 0.20t°:}^ and is also consistent with the quasar redshift, though 
with larger uncertainties. Neither fit provides evidence for extreme dust reddening as might be expected 
for a sub-mm-detected galaxy, and neither does a fit using solar metallicity spectral templates. Thus the 
identification of Q 0835+580 (117) with Q 0835+580 (SMI) is doubtful. However, such simple positional 
identifications of optical counterparts to sub-mm sources are not expected to be extremely robust given the 
14" FWHM of the SCUBA beam at 850/um. 

Figure lib shows that Q 0835+580 (458) is best fit as an instantaneous burst at Zp^ = 1.2t_Q'1 with 
E{B — V)= 0.55to;o5' though a constant star formation history at Zph = 1.35to!i5 with the same reddening 
is also acceptable. This galaxy's photometric redshift, which places it slightly foreground to the quasar, is 
consistent with it being slightly brighter than expected for a brightest cluster galaxy at the quasar redshift. 
Its high E{B — V) is consistent with it being the coimterpart of Q 0835+580 (SM2). This possible luminous 
sub-mm-detected starburst at (or nearly at) the quasar redshift clearly warrants spectroscopic foUowup. 



7.2. Q 1126+101 Field 

The field of Q 1126+101 (zq = 1.5173) has data in five filters (rzJHK) for most objects, and we can 
fit only objects with data (including upper limits) in all filters. There is no distinct spatially compact group 
of galaxies around Q 1126+101 as there is around Q 0835+580. However, there is a nearby clump of J — if 
selected EROs which may be galaxies background to the quasar or very dusty galaxies at the quasar redshift 
(discussed in §7.2.2), as well as a population of R — K selected EROs which may be old galaxies in an 
overdensity at the quasar redshift (discussed in §7.3). First, however, we discuss a possible sub-mm source 
identification. 



7.2.1. Q 1126+101 Field Possible Sub-mm Counterpart 

Q 1126+101 (384) was chosen for SED fitting since it is a possible counterpart of sub-mm source 
Q 1126+101 (SMI) (§5). The best fit (x^=2.78) is an instantaneous burst with Zph = l-50+g;5^ and E{B - 
V) = 0.50tojl (Figure 12); the constant SFR fit also requires a large E{B - V) = O-SSlfJ;??. These large 
E{B — V) values support the idea that Q 1126+101 (384) is the optical counterpart of the SCUBA source 
in this field. The galaxy has K = 17.7, making it ~0^4 brighter than the magnitude of a brightest cluster 
galaxy at the quasar redshift as estimated by HG98, but it is within the observed scatter in such objects' 
magnitudes at z ~ 1 (Thompson, Aftreth & Soifer 2000). As with Q 0835+580 (458), this possible luminous 
sub-mm-detected starburst at the quasar redshift clearly warrants spectroscopic foUowup. 



7.2.2. Q 1126+101 Field J-K selected EROs 

Four of the brightest J-K selected EROs in the Q 1126+101 field had data sufficient for SED fitting. 
Figure 13 shows the best fit of these, Q 1126+101 (447) and Q 1126+101 (424). Q 1126+101 (447) is best fit as 
a dusty instantaneous burst consistent with the quasar redshift {zph = I-QIIq'^, and E{B — V) = 0.75t'^''i^). 
Q 1126+101 (424) is best fit as a less dusty, constant SFR model background to the quasar at >99.9% 
confidence {zph = 3.00^0^1 and E{B -V) = 0.46lg j^). This agrees with our suggestion in §7.3 of HG98 
that it is an example of a dusty star-forming galaxy at ^; > 2.5, with the 4000 A break in J or beyond. 
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Thc other two J - K selected EROs we fit were Q 1126+101 (381), which has very red 2 - J > 3.4, 
and Q 1126+101 (425), the brightest J-band dropout in the field. In HG98 we suggested that these were 
examples of dusty star-forming galaxies at the quasar redshift Zq = 1.5173, which puts the 4000 A break 
between z and J. Figure 14 shows that both of these objects are best fit with highly reddened constant 
SFR models (though for Q 1126+101 (381) an instantaneous burst is almost as good) background to the 
quasar at almost exactly 90% confidence. Q 1126+101 (381) has z^h = l-95;o;35 and E{B -V) = Q.hhtl-fi- 
Q 1126+101 (425) has Zpu = 2.35;jj;^^ and E{B -V) = 0.55l^;o?. This Zpu for Q 1126+101 (425) is slightly 
lower than the value of Zph = 3.5 ± 0.5 reported in Hall et al. (1999) based on earlier modelling, but the 
E{B — V) is the same within the uncertainties. The fits for both objects have large xt values because they 
underpredict the flux in the r band, and the same is true to a lesser extent for Q 1126+101 (424). This 
may represent a limitation of our simple models with uniform E{B — V). The addition of a small amount of 
unreddened light from a young stellar population could bring the fits into agreement with the r-band data 
while not destroying their agreement with the observed fluxes at longer wavelengths (cf. HG98). 

As mentioned in HG98, the large E{B — V) values for these objects are consistent with their J — K 
colors being comparable to, or redder than, the prototypical dusty ERO HRIO, which requires dust reddening 
to fit its SED at its known z = 1.44 in any reasonable cosmology (Hu & Ridgway 1994; Graham & Dey 
1996; Cimatti et al. 1997). However, the prediction of HG98 that Q 1126+101 (381) and Q 1126+101 (425) 
would be at the quasar redshift seems erroneous based on their photometric redshifts, though the quasar 
redshift is excluded at only 90% confidence. The Zph of these objects are plausible in the sense that their 
K magnitudes would be fainter than those of powerful radio galaxies (typically the brightest galaxies at 
any redshift) at the same redshifts (van Breugel et al. 1998). Their being background to the quasar is also 
consistent with the lack of a strong concentration of excess galaxies around the quasar on the sky in this 
field (§2.3). We also fitted solar metallicity models to these four objects to see if their E{B — V) values or 
redshifts depended strongly on metallicity. The fits were of comparable but slightly worse xt average, 
and none of the best-fit redshifts or E{B — V) values changed enough to affect our conclusions. 

In summary, fits to the SEDs of these four J — K selected EROs all require strong reddening, with 
E{B — V) values of ^ 0.6+0.3. Two are better fit by constant SFR models, while two have nearly identical 
for either constant SFR or instantaneous burst models. The best-fit photometric redshifts are 1.9 < Zph < 3, 
and only Q 1126+101 (447) is consistent with being at the quasar redshift at 90% confidence. Small amounts 
of unreddened light from a young stellar population are needed to match the r-band data for all three of 
the Zph > Zq objects. All these conclusions hold for both solar and 20% solar metallicity templates. These 
J — K selected EROs may be representative of a population of dusty and star-forming galaxies missing from 
optically selected samples. 



7.3. Very Red Objects and Extremely Red Objects 

Photometric redshifts were also computed for all objects with r — K > 5.4 and K < 19.6 in both 
fields. Such objects could be old stellar populations at high redshift or dusty starbursting galaxies, both 
populations of great interest to the study of galaxy evolution. Objects with r — K > b.A but r — K < 6.322 
are not red enough to meet our definition of an ERO, but objects with r — K > 5.322 are called EROs by 
Cimatti et al. (1999) and Daddi et al. (2000), and objects with r — K > 5.622 are called EROs by Pozzetti 
& Mannucci (2000) and Martini (2001). To avoid confusion with our definition of EROs we refer to objects 
with 5.4<r-K < 6.322 as Very Red Objects (VROs; cf. Cohen et al. 1999). There are 12 such objects in 
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the Q 0835+580 field, one of which was not successfully fit . Excluding two objects with the colors of late M 
or L dwarf stars (cf. Hall 1998), there arc 12 VROs in the Q ff 26+101 field with sufficient photometric data 
for fitting, one of which was not successfully fit. Overall, in both fields, we successfully Hi 22 r — K > 5.4 
objects with K < 19.6.^^ These 22 objects include four from the spatially compact group discussed in §7.1.1 
and three of the four J-K selected EROs discussed in §7.2.2. One of these 22, Q 1126+101 (91), is best fit 
as a <0.2 Gyr old instantaneous burst with E{B — V)= 0.6 ± 0.2 at Zph = 5.0 ± 0.3. If this is correct, with 
K = 19.5 it would be one magnitude brighter than the z = 5.19 radio galaxy TN J0924-2201 (van Breugel 
et al. 1999), but consistent with the scatter seen in the K — z relation at 1 < ^; < 3. Thus this Zph is not 
impossible, though it is unlikely. In this section we discuss only the remaining 21 of 22 very or extremely 
red objects, which all have 1 < Zph < 2. 

Figure 15a plots E{B — V) vs. Zph for these 21 objects, and Figure 15b plots age vs. Zph- The 
distribution of E{B — V) values, and to a lesser extent the distribution of ages, are both skewed to slightly 
larger values than in Figure 8, with median E{B — 1/) ~ 0.4 instead of ~ 0.25 and median age ~ 0.7 Gyr 
instead of ~ 0.3 Gyr. Both trends are understandable since the objects in Figure 8 were chosen only for 
their proximity to the quasar, while those in Figure 15 were chosen on the basis of their red colors. Figure 15 
also shows that there are more objects consistent (at the 90% confidence level) with the quasar redshift in 
the Q 0835+580 field (open squares) than in the Q 1126+101 field (filled squares). This is expected, since 
many of the Q 1126+101 field objects are J — K selected EROs, and at least some of that population is not 
expected to be at the quasar redshift (§2.3 and §7.2.2). 

Following Pozzetti & Mannucci (2000), we plot these 21 objects on an r — if vs. J — K color-color 
diagram to study the separation of unreddened and reddened galaxies at 1 < Zph < 2 (see also Martini 2001). 
Figure 16 shows most objects' fits require dust reddening of E{B — V)> 0.3 (filled symbols, with symbol 
size scaling linearly with E{B — V)). The separation of reddened and unreddened galaxies (open symbols) 
is not as clean as predicted by the modelling of Pozzetti & Mannucci (2000) for objects at 1 < z < 2. 
Our modelling is very similar to that of Pozzetti & Mannucci (2000), but we fit a wider range of ages and 
E{B — V) values and use information from more filters than just RJK or UK. In particular, they do not 
consider reddened instantaneous burst models, which we find are the best fits for many of our objects. Larger 
samples plus spectroscopic redshifts and spectral types are clearly desirable to calibrate this classification 
system. For now, we note that the general trend of objects with larger E{B — V) being redder in J — is 
valid, but that this trend is not accurately predictive on an object by object basis. E{B — V) — Q objects 
may be an exception, since they seem to occupy a distinct region of the diagram with little contamination. 

We also note that our results may differ from those of Moriondo, Cimatti & Daddi (2000), who assembled 
HST images of a heterogeneous sample of VROs and EROs in both blank and candidate high-redshift cluster 
fields. They found that VROs and EROs with irregular morphologies, which are likely to include a higher 
percentage of star-forming galaxies than objects with regular morphologies, tend to have redder R — K or 
I — K colors. We may be able to compare these irregular EROs to objects in our sample which have large 
best-fit E{B — V) values, which also should include a higher percentage of star-forming galaxies than galaxies 
with low reddenings. As seen in Figure 16 and discussed above, we find that objects with large reddening 
do not tend to have redder r — K colors, although they do tend to have redder J — K colors. Our heavily 



I'^The two unsuccessfully fit objects are Q 0835+580 (618) and Q 1126+101 (50). Q 0835+580 (618) has Zj,h < 0.9 
at 90% confidence, but its age and E{B — V) are unconstrained since the SED can be fit either as an extremely red 
z < 0.5 galaxy and an unreddened z ~ 1.5 galaxy. We suspect Q 1126+101 (50) could not be fit because of its red 
r — z and blue z — J colors, which can probably be simultaneously fit only by a solar metallicity model at ~ 1. 
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reddened EROs and the irregular EROs of Moriondo, Cimatti & Daddi (2000) may bo different populations, 
of course, but it would be very interesting if they were in faet comparable. In that case this result would 
indicate a breakdown in the relationship between morphology and star formation properties at z > 1; for 
example, if we are seeing irregular galaxies with relatively old or quiescent stellar populations, and young or 
actively star forming galaxies with nonetheless regular morphologies. 

To test the sensitivity of our fits to metallicity, we refit 16 objects using solar metallicity models instead 
of 20% solar ones. These solar fits had comparable xt^ with a tendency toward worse fits on average. Best- fit 
redshifts and E{B — V) values were unaffected by the different metallicity, with no discrepancies at >90% 
confidence. Best-fit ages were found to be more sensitive to the metallicity, but even so, only in five cases 
where the solar metallicity fits had lower xt did the best-fit ages change (decrease) at >90% confidence. 
All five cases arc best-fit instantaneous bursts, and four are E{B — V") = objects. This is consistent with 
our E(B — V) values being essentially the same for solar and 20% solar metallicity fits; the majority of our 
objects require dust for both metallicities, while a minority do not require it for either metallicity. 

Finally, we attempted to fit four objects with z — J > 2.9, suggestive of a spectral break near l/xm, 
which have neither r — K > 5.4 nor J — K > 2.5. One object — Q 1126-1-101 (253) — was not successfully fit 
due to its detection only in J and K. The three successfully fit objects ar(^ denoted 'red z — J' in the Notes 
column of Table 4, which summarizes the fit results. Q 0835-1-580 (343) and (379) have fits similar to those 
of VROs and EROs. The r — K > 5.4 definition of a VRO is somewhat arbitrary, and these two red z — J 
objects have r — ii' ~ 5, so the similarity of their fits to those of VROs is not surprising. Q 1126-1-101 (269) 
also has a best fit similar to those of VROs or EROs, but it has a large xl value. We believe this is because it 
has only a moderately red color r — K = 4.4 which our simple stellar population models have trouble fitting 
simultaneously with this object's very red z — J > 3.7. A higher metallicity model with a slight contribution 
from an unreddened young population is probably needed to improve the Xv ^or this object. 

In summary, we have performed SED fits to a sample of red objects with a high a priori probability of 
having strong dust reddening or old stellar populations. Many of these objects' SED fits yield photometric 
redshifts consistent with the quasar redshifts. We find that significant dust is required to fit most of these 
objects, including about half of those which also require old ages. Our fits have E{B — V) values ranging from 
to 1.1, with median E{B — V) ^0.4. For comparison, Lin et al. (2001) find an average E{B — V) ~ 0.1 
from fitting the luminosity density evolution of z < 0.65 field galaxies from the CN0C2 redshift survey, 
Thompson, Weymann & Storrie-Lombardi (2000) find a distribution E{B — V) < 0.6, albeit peaked at 
E{B ~V) < 0.1, for 2; < 2 galaxies in the deep NICMOS image of the HDF-N, and Hammer et al. (2001) 
find E{B — T^) < 1, peaked at E{B — V) ^ 0.5, for a sample of luminous compact galaxies a± z ^ 0.65. 
Our E{B — V) values do not change significantly when we fit solar metallicity models instead of 20% solar 
models. However, the best-fit ages for the four objects with best-fit E{B — V) =0 do decrease significantly 
for solar metallicity fits. 

Thus the very red objects in our fields appear to consist of two populations. The minority population 
has no dust and is red due to age and/or high metallicity, while the majority population is red due to dust. 



^^This object is one of seven in Table 4 with xt ^ 10- Two others were discussed in §7.2.2; they also require a 
contribution from an unreddened young population to fit the observed SED. Of the remaining four, Q 0835-1-580 (307) 
probably requires both a metallicity > 0.2Zq and a young spectral component while Q 0835-1-580 (280) and (410) 
probably require only metallicities > Q.2Zq and Q 0835-1-580 (66) is definitely better fitted by a solar metallicity 
model. It is encouraging that all seven objects' high can be understood as due to limitations of our modelling. 
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albcit with old age contributing to the red colors in some cases. This result may differ from those of Stiavelli 
& Trcu (2000) and Moriondo, Cimatti & Daddi (2000), who find that the fraction of dusty starbursts among 
VROs and EROs is <30% and ~20— 50%, respectively, based on their morphologies as seen by HST. As 
discussed above, our heavily reddened objects and the irregular and exponential disk objects of Stiavelli & 
Treu (2000) and Moriondo, Cimatti & Daddi (2000) may be different populations, but it would be very 
interesting if this discrepancy was real and indicated a breakdown in the relationship between morphology 
and star formation properties at z > 1. 



8. Conclusions 

In previous work we have identified an excess population of predominantly red galaxies around a sample 

of 31 radio-loud quasars at 1 < z < 2. In this paper we have presented new multiwavclcngth data and analyses 
on the fields of four of these quasars at z ~ 1.54, obtained to build more detailed pictures of the environments 
of these quasars and the galaxies within them. Our conclusions are as follows. 

1. These fields have a surface density of extremely red objects (EROs, with R — K>Qi) 2.7 times higher 
than the general field. Assuming these EROs are passively evolved galaxies at the quasar redshifts, we find 
that they have characteristic luminosities of only ~ L* . Thus these r — K selected EROs are consistent with 
being drawn from the bright end of the luminosity function of early-type galaxies. The higher density of 
such EROs in RLQ fields can be easily understood if the RLQs are located in overdense regions. However, 
evolution in the dust content of these EROs could affect the luminosity estimates in up to ^40% of the 
population. Wc also show that only one of four RLQ fields has an excess oi J — K selected EROs with 
J — K > 2.5. The majority of these EROs arc therefore probably unrelated to the quasars. 

2. Wide-field J and Ks data show that the galaxy excess around Q 0835+580 extends to 140" at 2.35cr 
significance, with a richness of A'^o.5=27±ll (corresponding to Abell richness 2±1). The galaxy excess around 
Q 1126-1-101 extends to only 50", suggesting that this quasar is embedded in a small-scale overdensity even 
though the overall counts in the field are higher than the literature average, 

3. In three fields totaling 10.156 arcmin^ we present the deepest narrow- band rcdshiftcd Ha observations 
published to date. Wc detect five candidate galaxies at the quasar redshifts, a surface density 2.5 times higher 
than in the only existing random-field survey of similar depth. However, photometric SED fitting of one 
candidate suggests it is background to the quasar, and that [O iii] A5007 instead of Ha was detected in the 
narrow-band imaging. 

4. Submillimeter observations of three fields with SCUBA detect two of the quasars and two galaxies 

whose SEDs are best fit as highly reddened galaxies {E{B — V)^ 0.55) at the quasar redshifts. While many 
galaxies whose SEDs indicate the presence of considerable dust are not detected, the SCUBA limits are only 
sufficient to rule out the hypothesis that these galaxies are hyperluminous infrared galaxies. 

5. if -band adaptive optics imaging is used to estimate structural redshifts for two moderately red 
bulge-dominated galaxies in the Q 0835-1-580 field using the Kormendy relation between central surface 

brightness and half-light radius. Both objects have structural redshifts consistent with ellipticals foreground 
to the quasar at 2 < 0.2 or 1 < z < 1.35. We have calculated photometric redshifts for these two objects for 
comparison to the structural redshifts. One object is not well fit by any of the simple models we consider, 
whereas the other is only fit by a moderately reddened and young galaxy at Zph = 4.50 ± 0.15. It seems 
likely that our optical photometry for these objects is corrupted by scattered light from the nearby bright 
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star used for adaptive optics, resulting in inaccurate SED fits. Optical adaptive optics or HST imaging may 
be necessary to obtain accurate photometry for future study of these objects. 

6. Quantitative SED fits and the resultant photometric rcdshifts arc presented and discussed for numer- 
ous galaxies from specific populations of interest. Thirteen of eighteen objects in the spatially compact group 
around Q 0835+580 are consistent with being at the quasar redshift at >90% confidence. One of the thirteen 
is a candidate for a very old galaxy with no ongoing star formation; the others appear to have some ongoing 
star formation, as indicated by a best-fit constant star formation model or a young best-fit age. Fits to four 
J - K selected EROs in the Q 1126-MOl field ah require large reddenings of E{B -V) ~ 0.6 ± 0.3. One is 
consistent with the quasar redshift at 90% confidence, while the remaining three have best-fit photometric 
redshifts of 1.9 < Zph < 3. These objects may be indicative of a population of dusty, star- forming high- 
rcdshift galaxies potentially underrepresented in optically selected samples. Fits to 21 very or extremely red 
objects show that many of them have photometric redshifts consistent with the quasar redshifts. Significant 
dust is required to fit most of these objects, including about half of those which also require relatively old 
stellar populations. 

Overall, our observations support the hypothesis that radio-loud quasars at z > 1 can be found in 
galaxy overdensities. Among all but the very reddest galaxies in these overdensitics, ongoing or recent star 
formation with moderate amounts of dust seems to be common. A similar result has been obtained by 
Yamada et al. (1997) for the galaxy excess around the RLQ 1335.8+2834 sk, z = 1.1, and a variation in 
recent star formation histories has also been suggested for the large scale structures at z ~ 1.2 in the field 
of 3C 324 (Kajisawa et al. 2000). These results suggest that the z> 1 redshift range is at least approaching 
the redshift range in which the majority of early-type cluster galaxies are undergoing significant bursts 
of star formation. However, our finding that most very or extremely red objects require significant dust 
reddening may conflict with HST imaging results that show only '^20—50% of such objects have disklike or 
irregular morphologies suggestive of recent or ongoing star formation. This raises the intriguing possibility of 
a breakdown in the relationship between morphology and star formation properties at 2 > 1. Spectroscopy 
will be needed to confirm these speculations, of course, and is needed in any case to confirm or calibrate our 
SED fitting results. 
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This preprint was prepared with the AAS L^T^^X macros v5.0. 



Table 1. Observations 



Quascir 


Q 0835+580 


Q 1126+101 


Q 2149+212 


Q 2345+061 


Right Ascension (J2000) 


08:39:06.459 


11:29:14.19 


21:51:45.874 


23:48:31.836 


Declination (J2000) 


-1-57:54:17.12 


+09:51:59.6 


+21:30:13.51 


+06:24:59.25 


Redshift 


1.5358±0.0006 


1.5173+0.0010 


1.5385+0.0008 


1.5396+0.0012 


Wide-Field Observations 


Telescope+Instiument 


SO 2.3m+PISCES 


MDM 2.4m+TIFKAM 






J Date 


1999 Jan 6 


1999 AjH- 6 






J Exposure 


10560 


10560 






J 3<T Limit 


22.51 


22.37 






Ks Date 


1998 Sep 30, 1999 Apr 27-28 1999 Apr 7 






Ks Exposiare 


8980 


12360 






Ks 3(T Limit 


20.63 


21.03 






Narrow-band Ha Observations 


Telescope+Instrument 


IRTF+NSFCAM 


CFHT+REDEYE 


IRTF+NSFCAM ; CFHT+REDEYE 


Date 


1997 Mar 19, Oct 26-30, Nov 1 


1993 Sep 9 


1997 Oct 26-30, Nov 1 ; 1993 Sep 9 


Filter A 


1.6549-1. 6734^tm 


1.6631-1. 6693Atm 


1.6574-1.6759Mm ; 1.6631-1.6693/:im 


H-haxiA Exposure 


7920 




810 


5940 ; 570 


Ha Exposure 


34620 




4200 


37440 ; 7200 


Area 


0.926 




3.982 


0.697 ; 4.551 


f/fa •i'^ limit'' 


1.71 




0.78 


1.92 ; 0.87 


JCMT + SCUBA Observations 


Date 


1999 Mar 26 


1999 Mar 25-26 




1998 July 23 


Exposure 


7040 


15360 




10240 


TCSO 


0.02-0.08 


0.02-0.07 




0.055 


RMS Noise @ 850/im 


2.21 mjy/beam 


1.83 mJy/beam 




4.2 mJy/beam 




CFHT Pue'o 


+ KIR i?-band Adaptive Optics 


Observations 




Date 


1999 Dec 30-31 








Exposure 


10350 









Units of 10 ^"^ ergs cm ^ s ^ . 



Note. — Redshifts are from the detailed study of Tytler & Fan (1992), except for Q 1126+101 which was computed using their method. All observation 
dates are UT. All exposure times are in seconds. All areas are in arcmin^. 



Table 2. Candidate Ha Emitters 



ID H {H-Ha)AB Ha REW^ ind^ SFR^ zho!^ ID#^ Aa A5 



Q 0835+580 (Hal) 


20.13+0.17 


1.56+0.25 


316+79 


13.9+1.3 


14.7+1.4 


1.5358+0.0141 


398 


4.12 


22.31 


Q 2149+212 (Hal) 


20.21+0.18 


1.17+0.25 


50.2+18.9 


1.5+0.2 


1.6+0.2 


1.5389+0.0047 


171 


12.48 


-21.66 


Q 2345+061 (Hal) 


18.27+0.01 


0.54+0.07 


51.6+3.8 


4.6+0.2 


4.9+0.2 


1.5396+0.0141^ 


195 


-13.86 


28.06 


Q 2345+061 (Ha2) 


>20.63 (2cr) 


1-781^53 


1141- 


1 7+0.6 
■■-•'-0.4 


1.8+0.5 


1.5389+0.0047 


133 


8.96 


-1.08 


Q 2345+061 (Ha3) 


>20.63 {2a) 






J-^-0.4 


1 q+0.7 

-^•^-0.4 


1.5389+0.0047 




22.79 


21.47 



^Rest Equivalent Width in units of A, calculated from Eq. 5 of Bunker et al. (1995). 

^Units of 10~^^ ergs cm"^ s~^, calculated from Eq. 6 of Bunker et al. (1995) using the appropriate zeropoint. 

^Star Formation Rate in units of Mq yr~^ for Ho=75 km Mpc~^, 0^0=0.1, derived from the relation of Kennicutt (1983). 

^Possible redshift range determined from the bandpass of the filter in which the excess was detected. 

^ID number in our catalog of the field (available on request from the first author; cf. Table 4.) 

'^Excluding the range 1.5389+0.0047 (see text). 



Note. — 



Aa and A^ are offsets in arcseconds from the quasar in that field, with positive offsets to the North and East. 



Table 3. Probable and Possible Sub-mm Detections 



Son TCP 




Aa 


A6 


S/N 




CC 


Catalog; ID 


Note.s 


Q 2345+061 










2.8 


11.7 




157 




Q 0835+580 


(SMI) 


8.5 


-44.8 


4.6 


10.1 


0.37 


117 


tentative 


Q 0835+580 


(SM2) 


14.7 


36.1 


3.4 


7.5 


0.29 


458 


tentative 


Q 1126+101 




-2.0 


-0.7 


3.4 


6.2 


0.49 


350 




Q 1126+101 


(SMI) 


-14.7 


14.7 


3.0 


5.5 


0.46 


384 




Q 1126+101 


(SM2) 


-45.7 


-13.7 


4.4 


8.0 


0.71 






Q 112G+101 


(SM3) 


17. 1 


70.6 


G.7 


12.6 


0.69 




(outatiw: odgo 



Note. — Act and A6 arc offsets in arcseconds relative to the optical position of the quasar. CC is 
the correlation coefficient for each source. Catalog ID is the ID number of the possible optical/near-IR 
counterpart (s) in our catalogs of these fields (available on request from the first author; cf. Table 4). 



Table 4. SED Fitting Results 



IDAaA5 K r-K J-K z-J Model xt Zpu E{B-V) Age Notes 



Q 0835+580 Field Objects 
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1, 
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ERO, radio hotspot 
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.2 


-7 


.3 
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.13 


1 
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2 
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13 
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2, 
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.2 


9, 


,6 
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2, 
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14 
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1, 
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1, 
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1, 
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1, 
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3, 
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Q 0835-1-580 (Hal) 
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Table 4 — Continued 



IDAa AS K r-K J-K z-J Model xl %h E{B-V) Age 



Notes 



Q 1126+101 Field Objects 
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Note. — Catalog ID is the ID number of the possible optical/neax-IR counterpart (s) in our catalogs of these fields (avciilable on request from the first author). 
Aa and Ai5 are offsets in arcseconds from the quasar, with positive offsets to the North and East. Magnitudes are FOCAS total magnitudes; colors are FOCAS 

isopliotal colors. All lower limits to colors arc 3cr limits. The Model column gives the best-fit of the two star formation histories we considered: an instantaneous 
burst (Burst) or a constant star formation rate (Const), both with 20% solar metallicity. Ages are given in Gyr. In the Notes column, ERO denotes r — K > 6.322 
and VRO denotes 5A<r-K < 6.322 (§7.3), red 2 - J denotes an object with 2 - J > 2.9 but r - < 6.322 and J - K < 2.5 (§7.3), and SCO denotes a member 
of the spatially compact group around Q 0835+580 (§7.1.1). 
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Fig. 1. — The radial distributions of all K < 20.5 galaxies relative to the two quasars with new wide- 
field JK imaging are plotted as points with error bars from Gehrels (1986). The average K < 20.5 
counts and ztltj RMS scatter from the literature data compiled in Hall & Green (1998) are plotted 
as dashed and dotted lines, respectively, (a) Q 0835+580 field. There is a clear galaxy excess at 
<20". The overdensity on larger scales appears to extend out to '-^140" at 2.35a significance, (b) 
Q 1126+101 field. There is an excess on all scales relative to the literature, but only a weak excess 
at <50" around the quasar relative to the counts in this field. The spike at ^150" is due to a galaxy 
grouping on the sky which by chance is aligned tangentially to the quasar. 
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Fig. 2. — Color- magnitude diagrams from IRTF H- and narrow-band imaging. Dotted lines enclose 
the expected range of objects with narrow-band excesses of < 3a significance, (a) Q 0835+580 field, 
(b) Q 2345+061 field. 
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Fig. 3. Color-magnitude diagrams from CFHT H- and narrow-band imaging. Dotted lines enclose 
the expected range of objects with narrow-band excesses of < 3a significance, (a) Q 2345-f061 field. 
The two candidate 3cr Ha emitters have H>20.6 and H—Ha >1.7. The three bright objects outside 
the 3cr significance range range are galaxies larger than the 2^'0 radius photometric aperture, (b) 
Q 2149+212 field. The candidate Ha emitter has H ~ 20.2 and H - Ha 1.2. 



Fig. 4.— CFHT REDEYE narrow-band Ha images, each 128" x 128". North is up and east is 
left. Bright objects produce elongated crosstalk "images" at the same position in the other three 
quadrants, (a) Q 2345-1-061 field. The candidate Ha emitters are marked, including one selected 
as a candidate only in IRTF data, (b) Q 2149-1-212 field. The quasar is the brightest object in the 
SE quadrant. The candidate Ha emitter is marked. 



Fig. 5.— K image of the Q 0835+580 field (data from HGC98) overlaid with 850/j,m SCUBA contours 
(this paper, 7.04 ksec exposure, la RMS noise 2.21 mJy/beam). The sub-mm and optical images have been 
aligned using the default JCMT pointing. Light (green) contours are positive (2.5, 3.0, 3.5... tr) and dark 
(red) contours are negative (-2.5, -3.0, -3.5... a), with the outermost dark (red) contour marking the edge 
of the full SCUBA maps. The quasar is marked by Q, and sub-mm sources discussed in the text by SM#. 
Positive features marked by a letter (or not marked at all) do not correspond to peaks in the cross-correlation 
map and so are not considered real. 



Fig. 6.— K image of the Q 1126+101 field (data from HGC98) overlaid with 850/;im SCUBA contours (this 
paper, 15.36 ksec exposure, la RMS noise 1.83 mJy/beam). Feature A, although strong, was seen only on 
night 2 and is probably noise. See Figure 5 for further details. 
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Fig. 7. — Redshift estimates for adaptive optics targets Q 0835+580 (106) (solid line) and 
Q 0835+580 (112) (dotted line) using the average rest frame S-band surface brightness-effective 
radius relation for giant elliptical galaxies. Labelled tick marks show the positions for various 
assumed redshifts, including every Az=0.1 from z=l to z=1.5. Points are local ellipticals from 
Sandage & Perelmuter (1990). Both curves assume go=0.225. Assuming qo=0.0 (0.5) would make 
little difference at z<0.52 but would shift the vertical curves rightward (leftward) by ~0.1 in the 
log by z=2. 
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Fig. 8. — (a) E{B — V) vs. Zph for the best fit model of twelve objects with K <20 (large symbols) 
and six with K > 20 (small symbols) within 21'.'5 of Q 0835+580. Instantaneous bursts (shown as 
squares) are preferred over constant star formation histories (triangles) in 9 of 12 cases at < 20 
but only 2 of 6 cases at K > 20. The vertical dotted line shows the quasar redshift. The error bars 
are projected 90% confidence limits (see text). Objects with E{B — V)=0 have been arbitrarily 
offset by small amounts in E{B — V) to show the error bars, (b) Age vs. Zph for the same objects. 
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Fig. 9.— SED fitting results for ERO Q 0835+580 (339). (a) 20% solar metallicity fit. (b) solar 
metallicity fit. In each plot, the left panel shows the observed flux points and the best-fit SED 
points connected by straight line segments. In all panels, the two different star formation histories 
we modelled are denoted by lines of different colors: gray for an instantaneous burst and black for 
a constant star formation rate. The middle and right panels are two perpendicular projections of 
the three-dimensional space of E{B — V), z and age. In each of these two panels, the cross and 
the plus sign mark the best-fit values for instantaneous burst and constant star formation models, 
respectively, while the contours show the projection of the 90% confidence volume. 
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Fig. 10.— SED fitting results for adaptive optics targets, (a) Q 0835+580 (106). (b) Q 0835+580 
(112). See Figure 9 for details. 
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Fig. 11. — SED fitting results for candidate sub-mm source counterparts in the Q 0835+580 field, 
(a) Q 0835+580 (117). (b) Q 0835+580 (458). See Figure 9 for details. 
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Fig. 12. — SED fitting results for candidate sub-mm source counterpart Q 1126+101 (384). See 
Figure 9 for details. 
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Fig. 13.— SED fitting results for J - ii' selected EROs. (a) Q 1126+101 (447). (b) Q 1126+101 
(424). See Figure 9 for details. 
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Fig. 14.— SED fitting results for J - K selected EROs. (a) Q 1126+101 (381). (b) Q 1126+101 
(425). See Figure 9 for details. 
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Fig. 15. — (a) E{B — V) vs. Zph for all successfully fit objects with r — K > 5A and K < 19.6 in 

the fields of Q 0835+580 (open squares) and Q 1126+101 (filled squares). The dotted vertical lines 
show the quasar redshifts. The error bars are projected 90% confidence limits (see text), (b) Age 
vs. Zpfi for the same objects. 
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Fig. 16. — r — K vs. J — K for all successfully fit objects with r — if > 5.4 and K < 19.6 in the fields 
of Q 0835+580 and Q 1126+101 (see Table 4). Triangles denote objects for which instantaneous 
burst models are the best fits, and squares objects for which constant SFR models are best. Filled 
symbols denote E{B — V)> 0.3, and open symbols best-fit E(B — V)< 0.3. Note that all constant 
SFR models have best-fit E{B — V)> 0.3. The symbol sizes scale linearly with E{B — V), with the 
largest best-fit value being 1.1. Objects without symbols are all best fit as E(B — V) = bursts. 
Error bars extending off the plot indicate 3a lower limits on the r — K or J — K colors. Objects 
above the dashed line are EROs by our definition {R — K > 6, using r = R + 0.322). The dotted 
lines show the elliptical (upper left quadrant) and dusty starburst (upper right quadrant) regions 
defined by Pozzetti & Mannucci (2000), again using r = i? + 0.322. Our SED fitting results indicate 
that their criteria are not as robust as might be hoped. 
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